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We consider the effects of critical superconducting fluctuations on the scaling of the linear a.c. 
conductivity, oiuj), of a bulk superconductor slightly above Tc in zero applied magnetic field. The 
dynamic renormalization-group method is applied to the relaxational time-dependent Ginzburg- 
Landau model of superconductivity, with <j{uj) calculated via the Kubo formula to O(e^) in the e = 
4 — d expansion. The critical dynamics are governed by the relaxational XY-model renormalization- 
group fixed point. The scaling hypothesis a'{Lo) ^ ^^~'*^^5'(aj^^) proposed by Fisher, Fisher and 
Huse is explicitly verified, with the dynamic exponent z ~ 2.015, the value expected for the d — 3 
relaxational XY-model. The universal scaling function S{y) is computed and shown to deviate only 
slightly from its Gaussian form, calculated earlier. The present theory is compared with experimental 
measurements of the a.c. conductivity of YBCO near Tc, and the implications of this theory for such 
experiments is discussed. 

PACS number(s); 74.25.Fy, 74.40.-^k, 64.60.Ht 



I. INTRODUCTION 



The discovery of high-temperature superconductors 
has, for the first time, made it possible to experimen- 
tally probe the critical region of the zero-field normal- 
superconducting transition since fluctuation effects in 
these materials are enhanced by the short coherence 
length and the high transition temperature Tc. It is nat- 
ural then to ask: If scaling and universality exist in the 
critical region, to which universality class does the tran- 
sition belong? From observations of the effects of critical 
superconducting fluctuations on thermodynamic proper- 
ties, such as the penetration depth pU2], magnetic suscep- 
tibility p™, specific heat p,p| and thermal expansivity 
M a consensus is emerging that the zero-field normal- 
superconducting transition is in the .static universality 
class of the three-dimensional, complex order-parameter 
(3D XY) model. In contrast, the effect of critical fluctu- 
ations on transport properties, such as the conductivity, 
depends on the nature of the dynamics near Tc and is 
much less explored. 

In general, conductivity measurements on high- Tc su- 
perconductors show an enhanced response above Tc due 
to the presence of superconducting fluctuations. Out- 
side the critical region this enhancement can be ex- 
plained in terms of the Aslamazov-Larkin H theory of 
non-interacting, Gaussian fluctuations, and its exten- 
sions P,pl]|. In these theories the dynamic exponent z 
associated with the growth of the characteristic order- 
parameter time-scale near Tc appears in the conductiv- 
ity and takes the value z — 2. By examining the de- 
viation of z from 2 inside the critical region through 
linear d.c. P,^,pd]-p^ , non-hnear d.c. |15-1^ and hnear 



a.c. pq | conductivity measurements, the dynamic univer- 
sality class can, in principle, be determined. Currently, 
however, there is much variation in the measured val- 
ues for z and the dynamic universality class of the zero- 
field normal-superconducting transition remains uncer- 
tain. Unlike d.c. measurements, measurements of the a.c. 
conductivity p3 can test the scaling of the conductivity, 
(t(lj), over a wide range of frequencies, uj, thereby pro- 
viding a stringent test of theory. In the experiments of 
Ref. [Q the a.c. conductivity exhibits a scaling collapse 
which deviates slightly from the Gaussian theory. How- 
ever, the Gaussian theory is known to break down in the 
critical region. Thus, to sharpen the comparison between 
experiment and theory, we go beyond the Gaussian de- 
scription of fluctuations in this paper and calculate the 
scaling behaviour of the a.c. conductivity in the critical 
region of strong, interacting fluctuations. 

Fisher, Fisher and Huse (FFH) ||l^ have argued that 
near a second-order phase transition, if dynamic scaling 
holds, the a.c. fluctuation conductivity should scale as 
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where the correlation length for fluctuations in the su- 
perconducting order-parameter at temperature T is ^ ~ 
\T — Tcl"'^ with the static exponent j^, d is the spatial 
dimensionality, z is the dynamic exponent and S{y) — 
S'{y) + iS"{y) is a universal, complex function of the 
scaled frequency y ~ f^C^, with real and imaginary parts 
S' and S" , respectively. Outside the critical region, and 
in the d.c. limit, Eq. (1.1) reduces to the Aslamazov- 
Larkin theory. Since the conductivity is causal, and 



also finite for non-zero frequencies, Eq. (1.1) leads to the 
power-law behaviour at Tc 



a{u;) ^ (-ic^)-(2-'^+-)/-, 



(1.2) 



reflecting the absence of a characteristic time-scale at 
criticality. At Tc the phase 



{lo) = tan 
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(1.3) 



of the conductivity is independent of frequency, with the 
value Hol 



TT f2-d+z 
2 



(1.4) 



Equations (P_j) and ( |l.4[ ) allow one to determine the dy- 
namic exponent z independently of the static exponent v 
through a measurement of the a.c. conductivity at criti- 
cality. However, to go beyond these two results and cal- 
culate the entire universal scaling function S{y) requires 
knowledge of the renormalization-group fixed-point that 
determines the universality class for the dynamics near 

The time-dependent Ginzburg-Landau (TDGL) model 
of superconductivity provides an appropriate framework 
in which to study dynamic critical behaviour in this sys- 
tem pO| , pi| . Since this is the first detailed study of the 
dynamics in the critical region of the superconductor, 
and given the uncertainty as to which dynamic universal- 
ity class describes the transition, we consider here only 
the simplest, relaxational, dynamics for fiuctuations in 
the superconducting order-parameter — model A in the 
Hohenberg and Halperin classification |2^j2^. Previous 
studies of this model have implemented the Gaussian ap- 
proximation, where quartic interactions among fluctua- 
tions in the Ginzburg-Landau free-energy are neglected 
P,10l^In this approximation, the conductivity scales as 



Eq. (O) with V = 1/2 and z = 2, the exponents for the 
Gaussian flxed-point, and the scaling function 5(w^^) has 
been explicitly calculated. 
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FIG. 1. Com par ison of the a.c. conductivity scaling func- 
tion S(y), Eq. (L.7), for the relaxational 3D XY critical the- 
ory (solid curve) with the scaling function, Eq. (3.7), for the 
Gaussian theory (dashed curve). To facilitate later compar- 
ison with experiments [[L8|, the magnitude of S{y) is plotted 
against the scaled frequency t/ on a log — log scale. 

In the critical region the Gaussian approximation 
breaks down since the quartic interactions become im- 
portant, producing the critical flxed-point for the relax- 
ational XY- model 123024]. In the e — A — d expansion, the 
exponents for this flxed-point are well known pO| . An ex- 
trapolation of the 0{e^) results to three dimensions gives 
a correlation-length exponent oi v ~ 2/3 and a corre- 
lation function exponent of 77 « 0.02. For relaxational 
dynamics the dynamic exponent z is, to O(e^) p4| : 



with 



cq 



6 In 4/3-1, 



(1.5) 



(1.6) 



giving z « 2.015 in three dimensions. 

In the critical region, and near four dimensions, we 
verify that the a.c. fluctuation conductivity satisfies the 
FFH scaling hypothesis (El) for the relaxational XY- 
model fixed-point. We compute the complex scaling form 
S{y) appearing in Eq. (1.1) to O(e^), with the result 



S{y) 



2z' 



1 



{d-2 + z){d~2)y^ 



d-2 + z 



^y-{l-^yy''-'+^y^ 



(1.7) 



where y ~ lu^^ and z is given by Eq. (1.5) with Eq. (1.6). 
This is the main result of this paper, and is the prod- 
uct of a much more involved analysis than that used to 
determine the exponent z. Sections p4[-|V]| provide the de- 
tails of the calculation. The result (1.7) has the scaling 
behaviour stated in Eq. (|1.2|). Since, to O(e^), z for the 
critical theory in three dimensions is only slightly dif- 



ferent than two, the scaling function S{y) for the critical 
theory is very close to t he G aussian result calculated ear- 
lier (see Fig. 1). In Sec. VI] we compare the experimental 
a.c. conductivity data of Booth et al. ||l^ to the critical 
theory, extrapolated to three dimensions, and comment 
in Sec. VIII on the implications of this work for such 



measurements. 



II. FORMALISM 

A. The time-dependent Ginzburg-Landau model of 
superconductivity 

We describe the critical dynamics of a superconductor 
with a complex order-parameter ijj using the relaxational 
time-dependent Ginzburg-Landau model 
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with the Ginzburg-Landau free-energy 



F 



'fr (iV^I^-frolVl' + ylV'l') 



(2.1) 



(2.2) 



In Eq. (2.1) Fq is the bare order-parameter relaxation 
rate. Both To and the bare coefficient uo, which appears 
in the free-energy (2.2), can be considered temperature- 
independent near the transition; however ro ~ T — Tco 
changes sign at the mean-field transition temperature 
Tco, becoming negative for temperatures below T^q- We 
choose units so that ti = kgTc — 1 and m = 1/2, where 
m is the mass of a Cooper pair. The superconductor 
is assumed to be isotropic. The complex noise field C 



in Eq. (2.1) is taken to have zero mean and correlations 
described by 



(C(r,i)C*(r',0>=2ro<5(r-r')<5(t-0, 



(2.3) 



where the brackets (• • •) denote an average over the noise 
distribution, assumed to be Gaussian. The factor 2ro 



in Eq. (2.3) follows from the fluctuation-dissipation the- 
orem and ensures that the system relaxes to the proper 
equilibrium distribution. 

We will work in the symmetric phase, T > T^, with 
zero applied magnetic field and consider order-parameter 
fluctuations about a mean of zero. Fluctuations of the 
vector potential are neglected |£^. Since we will use the 
Kubo formula to calculate the linear conductivity from 
the system in zero electric field, an electric field is not 
included in Eqs. (2.1-2.2). In the classification scheme 



of Hohenberg and Halperin [|T|, Eqs. (2.1-pT3) consti- 
tute model A dynamics for a two-component (complex) 
order-parameter. Thus our model is in the dynamic uni- 
versality class of the relaxational XY-model |^,^ . 

Since the Ginzburg-Landau theory is coarse-grained, it 
contains an ultra-violet (UV) cutoff, A (corresponding, 
for example, to the lattice constant) |g^. This cutoff is 
manifest in the definition of the Fourier transform of the 
order-parameter , 



(2.4) 



iP{r,t)= / ij{k,u;)e''^-^-^-'. 

J kuj 

For convenience, we employ the short-forms 



' d'^k 
doj 



(2^) 



(2.5) 
(2.6) 



for the wavevector and frequency integrals, with the 
wavevector integral restricted to |k| < A. The existence 
of the cutoff will be crucial when we interpret the results 
of the theory. 

The order-parameter correlation function and the re- 
sponse function are central in what follows. The order- 
parameter correlation function, C(k, w), is defined as 



C(k,c.) = (^(k,c.)^*(k,c.)). 



By adding a source term. 



Fh^ 



d'^r {h*iP + hiP*), 



(2.7) 



(2.8) 



to the free-energy (2.2) we can define the (linear) re- 
sponse function, G(k, w), as 



G{k,uj) 



6{lP{k,L0)) 



Sh{k,Lu) 



(2.9) 



/i=0 



This measures the response of the order-parameter to the 
source h. Near equilibrium, the correlation and response 
functions are related though the fluctuation-dissipation 
relation 12^, 



C{k,L0) 



Im G(k,w). 



(2.10) 



B. The Kubo formula for the conductivity 

The linear a.c. conductivity, <t{uj), for an isotropic ma- 
terial can be deflned in terms of the current response, J 
(which includes normal and supercurrent contributions), 
to an inflnitesimal applied electric fleld, E, through 



3{uj) = cr(cj) E,{lu) 



(2.11) 



Since the quantities in Eq. (2.11) are evaluated at zero 
wavevector we suppress their wavevector dependence. 
The conductivity is complex and has a real dissipative 
response, a' , and an imaginary reactive response, a": 



a{uj) = a'{uj) +ia"{uj). 



(2.12) 



In linear response, the conductivity is related to a cur- 
rent correlation function via the Kubo formula [ p8[ . Near 
Tc strong superconducting fluctuations give a singular 
contribution to the conductivity which dominates the 
non-singular contribution due to normal electrons. Thus 
we may use the Kubo formula to calculate the real part 
of the conductivity due to superconducting fluctuations 



from the supercurrent correlation function, evaluated at 
E = 0§: 

a'M = ^(J.H.J.(-^))|E=o. (2.13) 

The supercurrent, J^ is 

J,(r, t) = -zeo(V* VV- - ^Vt/;*), (2.14) 

where cq is the bare charge of a Cooper pair. The imagi- 
nary part of the conductivity can be obtained by applying 
the the Kramers-Kronig relations pq] to Eq. ( 2.13| ) . 

The average in Eq. ( 2.13) is a four-point order- 
parameter average since Js ( p.l4 ) is quadratic in ip. 
Quite generally, this four-point average can be written 
as the sum of a "disconnected" product, ct^^-*, of two 
two-point averages, and a "connected" four point-average 



(4). 



a 



with 



cr'{uj)=cr^^>{uj)+a^'^>{uj) 



(2.15) 



2ei^ 
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a^-'\uo)^^ [ fc2 C(ki,c^i)C(ki,c^i-^w) (2.16) 



and 



„2 /.A 



^(4)(^) = ^ /■ ki .k2 C(4)(ki,..i,k2,u;2: 



W 



(2.17) 

where the exact two-point order-parameter correlation 
function, C(k, w), is defined in Eq. (2.7) and 



C^'')(ki,wi,k2,W2;w) = 

(V'(ki,CJi)V'*(ki,LJi -tj)V'(k2,W2)V'*(k2,t^2 +w))c 

(2.18) 

is the connected four-point order-parameter correlation 
function. 



C. Iterative dynamic perturbation theory 



The order-parameter averages (2.7) and ( 2.18| ) that ap- 
pear in Eqs. ( 2.16| ) and ( 2.17| ) can be expanded as a 
perturbation series in the bare non-linear coupling uq 
appearing in Eq. (p^). Dynamic perturbation theory 
for the time-dependent Ginzburg Landau equation (2.1) 
can be implemented either by using a Martin-Siggia-Rose 
p9| , p0[ field-theoretical formalism, or by a direct iteration 
of the equation of motion i2Q,pl| . The iterative approach 
involves less formal machinery and will be used here. 
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FIG. 2. The diagrammatic representation of the equa- 
tion of motion (2.1!^). Wiggly lines correspond to the order- 
parameter i/j (a starred wiggly line is ill*). The dotted line 
represents the Gaussian field t/jq. The Gaussian response func- 
tion Go (2.21) is shown as a line with an arrow. The vertex, 
where the response function meets three wiggly lines contains 
a factor — uo, as well as V (2.22), which conserves wavevector 
and frequency at the vertex. Iteration corresponds to replac- 
ing the wiggly lines on the right-hand side with either the 
first or second term on the right-hand side. In this way, one 
generates a series in uq. 



The equation of motion (2.1) can be explicitly written 
in Fourier space as 

V'(k,w) = '0o(k,u;) -Uo G'o(k,a;) 

X / Vx l/'(ki,Wi)V'*(k2,W2)-(/'(k3,t^3), 

(2.19) 



where 



V'o(k,c^)-— Go(k,c^)C(k,c^), (2.20) 

t 



(2.21) 



Go(k,c^)= ( -—-^ro + fc^ 
to 



and 

V = {2tiY+H{\^ - ki -I- k2 - k3)(5(w -0Ji+UJ2~ ^3). 

(2.22) 

The Gaussian theory neglects the non-linear interac- 
tion (mq — 0). In this case Eq. (2.19) reduces to -0 = -00, 
and the order-parameter is a Gaussian field by virtue of 
Eq. ( p. 20 ) and the fact that C is Gaussian. The order- 
parameter correlation function (pj) can be evaluated us- 
ing Eq. (O) and is 



Co(k,o.) 



(0o(k,c.)0*(k,c.)) 

2ro 

c^2 + r2(ro + fc2)2- 



(2.23) 



If a term coming from Fh (2.8) is included in the equa- 



tion of motion ( 2.19|) , it is straightforward to show that 



Go (2.21) is the Gaussian response function. A glance at 
Eqs. (2.21) and (2.23) shows that the Gaussian theory 
satisfies the fluctuation-dissipation relation, Eq. (2.1C). 

Since -0 appears in the integral on the right-hand side 
of Eq. (2.19), this equation can be iterated to produce an 



expansion for ■0 in powers of the bare coupling constant 
Mo- Averages containing tjj are then expressed as sums of 
higher-point Gaussian averages over jpQ, which break up 
into products of Cq 's. To keep track of the algebra, it is 
helpful to use the graphical representation of Eq. ( 2.19| ) 
shown in Fig. 2. In the graphical context iteration corre- 
sponds to "putting branches on the tree" and averaging 
corresponds to joining two conjugate dashed lines (V'o) to 
form a correlation function Cq. By examining all possi- 
bilities for joining for a given average, a series of graphs is 
generated with the proper symmetry factors. In dynam- 
ical perturbation theory there are two propagators: the 
response function Gq, denoted by an arrow, and the cor- 
relation function Cq, denoted by a line with a circle on it. 
Wavevector and frequency are assigned to these lines on 
the basis of conservation of wavevector and frequency at 
the graph vertices, given by V in Eq. ( 2.22| ). Wavevectors 
and frequencies flowing around loops are integrated over. 
More details of the graph rules can be found in |^,|l| . 
An example of this procedure is the self-energy diagram, 
Fig. 3, and the corresponding algebraic expression (4.5). 



where the renormalization constant Z^ comes from 
"wavef unction" renormalization (a rescaling of "0) and, 
in the minimal subtraction scheme, is given by [ p2| , p3[ 

Z^ = l~-u^ + 0{u^). (2.28) 

The renormalized "mass" r is defined as 



rEEG^i(0,0), 



(2.29) 



which, using Eqs. (|2.2lD , ( |2.26D and ( p.27| ), is related to 
the bare mass tq by 



= Z^[ro-S(0,0)]. 



(2.30) 



Near Tc the physical response function at zero wavenum- 
ber and frequency behaves as Gu{0, 0) = ^'^~^k~'^, where 
r] is the usual correlation function exponent and ^ is the 
order-parameter correlation length which diverges as 



^ ~ |T - T, 



(2.31) 



D. Renormalization of the theory and the XY 
fixed-point 



with the critical exponent v. Thus, from (2.29), we have 
r = ^-2+"^". (2.32) 



It is computationally convenient to dimensionally regu- 
larize the theory and renormalize via minimal subtraction 
|^,Q. This will produce a sensible e — A — d expansion. 
To be more concrete, we define the renormalized "cou- 
pling constant," u, in terms of the bare coupling constant, 
Mo, by 



Z,, 



uo, 



(2.24) 



and define the dimensionless, renormalized coupling con- 
stant, u as 



Sd 



2(27r)'^ 



UK 



(2.25) 



where k is an arbitrary wavevector scale and Sd is the 
surface area of the unit sphere in d dimensions. The 
renormalization constant Zu — 1 + 0{u) [Q. Since only 
u^ will appear in the conductivity, and we neglect terms 

of 0{u'^) and higher, we may approximate Zu ~ I- 

Renormalization of the bare response function ( |2.9| ) 
provides the remaining renormalization constants. The 
bare inverse response function including self-energy cor- 
rections, E, may be written 



G-^{k,uj) = Go i(k,cj) - E(k,cj). 



(2.26) 



The renormahzed inverse response function Gj^^(k, w) 
may be expressed in terms of the bare quantity ( 2.26| ) 
by 



G],\k,uj) = Z^G-\k,uj), 



(2.27) 



Since we are neglecting magnetic fluctuations and work- 
ing at the "uncharged" fixed point, the renormalized 
charge, e, is simply the bare charge: e — cq. Finally, 
the bare relaxation rate Fq, appearing in the dynamic 
response function ( ^.26 ) is related to the renormalized 
relaxation rate F by 



1 1 



(2.33) 



where, from minimal subtraction, the renormalization 
constant Zr for this relaxational model is |3flJ33] 



e 



(2.34) 



The constant c is given by Eq. (|1.6|). 

Near Tc, as one probes the long- wavelength physics, 
the coupling u flows towards the fixed point value u* de- 
termined by the IR-stable zeros of the renormalization- 
group beta function P{u*) = [Q. This mechanism 
is responsible for universality. To leading order in the 
e-expansion, u* is [p3,p3| 



10 



+ Oie'). 



(2.35) 



This is the Wilson-Fisher |23| fixed-point for the XY- 
model. The correlation function exponent 77 is related to 



Z^, (2.28) and has the following expansion in u*: 

r; = 2(r)2 + 0((r)3). (2.36) 



The result v sa 2/3 quoted in the Introduction, which also 
appears in ( ^.31 ) , is an extrapolation of the e-expansion 
result to three dimensions. Finally, the dynamic expo- 



nent z is related to Zr (2.34) for the relaxational dynam- 
ics, and given hy z — 2 + cr] with c = 6 In 4/3 — 1 and rj 



given by Eq. (2.36). Thus, by reorganizing the theory as 
an expansion in e and using the fixed point value u* for 
the coupling, the IR divergences near criticality can be 
sensibly treated and lead to corrections to the Gaussian 
exponents. 

Even after we renormalize the conductivity as de- 
scribed above, some poles in e will remain. These poles 
are due to UV-divergences in the theory for the conduc- 
tivity that appear even at the Gaussian level and have 
nothing to do with the critical behaviour. These poles 
must be eliminated by adding a constant to the conduc- 
tivity, as will be discussed in Sec. 



VI 



III. THE CONDUCTIVITY IN THE GAUSSIAN 
APPROXIMATION 

We now review earlier work on the a.c. conductivity in- 
volving non- interacting, Gaussian fluctuations P,pi3|, and 
set uq = in Eq. (2.2). In the Gaussian approximation 
the connected piece of the conductivity, Eq. (2.17), is 
zero. Thus, from Eqs. (2.15) and (2.16) one has 



2p2 /-A 
a'ico) = ^ kl Co(ki 

" Jfciwi 



wi)Co(ki,cji+a;), (3.1) 



where Cq is given by (2.23). The calculation of the inte- 
gral in Eq. (3.1) involves a contour integration over the 
frequency variable, and then a straightforward evaluation 
of the remaining wavevector integral, with the cutoff A 
set to infinity. The complex conductivity takes the form 
~T0l: 



e2 _ ^4-'* 
"^") = 2f^^ 43d ^«(2^")' 



(3.2) 



where 



Sd 



{2'kY 
and the scaled frequency j/o is 



r(d/2)r(3-rf/2) 



(3.3) 



yo 



2ro' 



The Gaussian order-parameter correlation length ^o is 
defined as 



C — -1/2 
f = r-Q , 



(3.5) 



thus ^0 "^^ |T — T'col^^/'^ and v = 1/2 in the Gaussian 
theory. The real part of the scaling form Sq is computed 
from Eq. (3.1) to be 



S'g^Vo) = 



d{d-2)yl 



1-(1 



vl 



)'*/4cos(|tan-iyo 



(3.6) 



The imaginary part of the conductivity is obtained from 
Eq. ( p.6| ) using the Kramers-Kronig relations. The result 
for the complex scaling form is then 



Sg{vo) 



1 



d{d-2)yl 



1 - ^ *yo - (1 - iVoY'^ 



(3.7) 



The Gaussian result, Eq. ( p^ w ith the definition (pl|), 
satisfies the FFH hypothesis (|l.l|) with z = 2. 

We note two properties of these results that will be 
important later. The first is that Eq. (3.2) has a factor 
of e = 4 — d in the denominator. This is a consequence 
of setting the cutoff A to infinity, and indicates that even 
the Gaussian theory is sensitive to the cutoff in four di- 
mensions. The second property is that S'q {'L&) has the 
e-expansion 



S'civ^) 



1 



Y.'' S^y,) 



i=l 



The coefficient of e in Eq. (|3j) 
3 1 



Siivo) 



2[(l-2/o)ln(l 



y'o) 



(3.8) 

4yotan~iyo], 
(3.9) 



is interesting because it appears later in both the discon- 
nected and the connected pieces of the conductivity. 



IV. DISCONNECTED PIECE OF THE 
CONDUCTIVITY 

To go beyond the Gaussian theory requires the calcula- 
tion of both the full two-point correlation function (^.7|), 
including self-energy corrections, and the four-point aver- 
which appear in the conductivity through Eqs. 



age (p.l 

( ^.16| ) and ( 2.17 ) . The calculations must be performed to 
0{u^), where the first corrections to the Gaussian result 
z = 2 occur. In this section we examine the disconnected 
piece of the conductivity ( 2.1(j ) . The next section tackles 
the connected piece. 

We first dimensionally regularize and renormalize the 
(3.4) theory as outlined in Sec. II D. From Eq. ( [2.16 ) 



the 



disconnected contribution to the conductivity is then 



r(2) 



(C) 



2e2 
d 



klC{ki,LUi)Ciki,u;i+io), (4.1) 



kiuJi 



where C is the full correlation function (2.7), includ- 
ing self-energy corrections. We will calculate the re- 
sponse function G (2.9) to 0{u^) and use the fluctuation- 
dissipation relation (2.10) to get C. With the definition 



( |2.30 ) of the renormalized mass, r, the inverse response 
function (2.26) may be written as 

G-\k,u;) = G^\k,u;) - [S(k,c.) - £(0,0)], (4.2) 

where now 



To ^ r/Zr 



(4.3) 
(4.4) 



in Go ( |2^ ) and Co (g^). To 0{u^} only the "Saturn" 
diagram Ss(k, lu) shown in Fig. 3 contributes to Eq. ( |4.2| ) 
since, to this order, it is the only piece of the self-energy 
that is wavevector- and frequency-dependent. Applying 
the rules outlined in Sec. uTQ to Fig. 3 gives 



E,(k,w) = 6u^ I Co(k2,U2)Coik3,LU3) 



I k2UJ2k3U3 

xG'o(k- k2 - k3,a; - UJ2- wj,). 

(4.5) 




FIG. 3. The Saturn diagram, Es, for the self-energy 
consists of two loops formed by two correlation functions 
Co (lines with circles) and one response function Go (line 
with an arrow). Wavevector and frequency flow through 



the diagram in accordance with the discussion in Sec. II C 



The correlation function C is then obtained from (2.1C) 
and (O): 



C(k,c^) = Co(k,c^) 



Im{G2(k,c^)p,(k,w)- 2,(0,0)]} + 0(^3). 

(4.6) 



Thus the disconnected piece of the conductivity (4.1) can 
be expressed in terms of the integrals 

h[Lo)^'^[ A:?Go(ki,c^i)Go(ki,c^i+c^) (4.7) 

and 

4e2 r 1 

-^2(w) = -rlm / /c2 Go(ki,cji) — Gg(ki,tJi+cj) 

x[E,(ki,c^i+co)- 2,(0,0)], (4.8) 
by writing 

a^^\uj) ^h{u) + 2l2{oj)+0{u^). (4.9) 

Each integral is dealt with separately below. 



A. The integral Ji 



The only differences between /i (4.7) and the starting 
point (3.1) of the Gaussian calculation are the substitu- 
tions: Eqs. (E^), (4.4) and eo -^ e. Transcribing the real 
part of the Gaussian result (3.2) gives 



h{u,)^^aK-^^zfZrS'a{y) 



with S'q given in Eq. (i.i) and 

uZ^Zy 



y = 



(4.10) 



(4.11) 



2Fr 

The dimensionless measure of the nearness to the transi- 
tion is 

x^^, (4.12) 

K 

where the arbitrary wavevector scale k was introduced 
earlier in Eq. (2.25). From the expression (2.32) for r we 
have 



m 



-l+v/2 



(4.13) 



The function SqIj)) can be expressed in terms of the 
scaled frequency y. 



y 



CJ^^K 



z^z-2 



2r 



(4.14) 



with z given by Eq. (1.5), by the expansion 

S'aiy) = S'aiy) + dyS'aiy) (y - y) + \dlS'a{y) [y ~ yf 

^\dlS'a{y)[y-yf + ---, (4.15) 

where d y ind icates a derivative with respect to y. The 
results (^.281) for Z^, ( ^.32[ ) for r and (|2.34| ) for Zy are 
used to obtain the following relation between y and y: 

y - y = y(^c + I) (tiIiix u^ ) +0{u'^e), (4.16) 



where c is give n by Eq. {^^^). Using equat ion (4.16 ), the 
expansion (3^) of S'q, and the fact that rj ( ^.3(j ) is O(e^), 
we write (4.15) as 



S'aiy) = 1 + eMy) + ^S2{y) + e^SM 

-fe(c+l) iri\nx~-uA [ydySi{y) + eydyS2{y)] 
+0(il2e2,e4). (4.17) 



We now use the expansions (2.28) of Z^ an d (2.34) of 
Zp together with Eq. (4.17) to write /i (4.10) as a series 
in w, with coefficients expanded in powers of e. Terms 
of 0{v?e, €^) and higher are neglected (since the fixed- 
point value u* (2.35) is 0(e) we are effectively working 
to 0(6^)). The result, written in a form that will be 
convenient for later analysis, is 



h{u;) 



2r 



1-2^' 



'^ -2 I 1 I '^-2i ''-2c (" ^ C-|-l_2 



;■« H 1 — M Inx M S'i(y) — 



6 — CU 

u^ ydySi{y) - ln.x H — (In x) 



6 



{lnxf + 



1 — (e — CU ) In a; + — (In a;) 



Si{y) + {e~ cu^){l - elnx) S2{y) + e^S^iy) 
+ {c+l){'n + u^)ydySi{y) Inx - {c+l)u^ ydyS^iy) + 0{uh,e^)}. (4.18) 



B. The integral I2 

The calculation of I2, Eq. ( |4.§| ), is in volved so we only outline it here. The first step is to re-scale the internal 
wavevectors and frequencies in Eq. ([4.8[) by 



LOi -^ Tr LOt, 
where i = 1, 2, 3 (remember that E^ contains an integral over \i.2'-^2^3^3) , and write 



(4.19) 
(4.20) 



(4.21) 



The dimensionless integral in Eq. ( [I.2l| ) 



/2(y) = 16Im / klCo{-ki,LUi) ^—-Gl{kuiUi + 2y)CoCk2,u;2)Co(k3,iU3) 

X [Go(ki - ka - kg, 2j/ -H t^i - W2 - W3) - Go(-k2 - kg, -W2 - t^a)], (4.22) 
is written in terms of the dimensionless functions 



and 



Go(k,c.) = 



Go(k,c^) = 



1 



l + fc2 



u;2 + (l + A:2) 



2^2 ■ 



(4.23) 



(4.24) 



Since I2 is already Oiv?') we have simply replaced all bare coefficients in Eq. (4.21) by rcnormalized ones, and used 
the scaled frequency y from Eq. (4.14). 

The second step is to evaluate the three frequency integrals in Eq. ( 4.22| ) by contour integration. The calculation 
is straightforward and yields 



/2(y)= Re[i^{y) + il{v) + F^{y)], 



with 



H^iy)^ / dv{l~v) ki 

Jo Jkik2k3 

pi 
Tbt„,\ o / J„. t^ „.\ I ;„2 



1 



a2as.{a2 + ay, + 04) (as + 2iy) 
1 1 



2 1 

a\ (ai + iyvY 

I',{y) ^Z I dv{l~v) I ki / 

Jo Jkik2k3 [ai + iyv)'^a2a3 

i^(y) = -3iy [ dv {1 - v) [ kf — 

Jo Jktk2k3 (fli +2yw)''a2a3 (02 + 03 +04) (ag + 2iy) 



02 + as + a4 02 + 03 + 04 
1 



where, for convenience, we define 



Oi = 1 + kf, 

04 = l + (ki+k2+k3)2, 

a4 = l + (k2+k3)2, 

05 = Oi + 02 + 03 + 04. 



i = 1, 2, 3 



(4.25) 

(4.26) 
(4.27) 
(4.28) 



(4.29) 
(4.30) 
(4.31) 
(4.32) 



Note that we have used the Feynman formula 



1 r{ai+a2) f\ ,. 

"sr-57 = -FT — vFT — V / ^^(1 



ai — 1 Q2— 1 



v)"^^^v 



[{1~V)C1 +VC2]"' 



+02 



(4.33) 



with the Feynman parameter v to group and simpUfy terms in Eqs. (4.26-4.28). 

The final step is to evaluate the wavevector integrals in Eqs. (4.26-4.28) using ( 4.33| ) and e-expand the resulting 
integrals over Feynman parameters to 0(e°). An example of this procedure appears in Appendix H. The results for 
Eqs. ([4.26[|4.28|) are 






1,4 1 , , ,, 4 0.003 _„, , ^, ' 
7T3 In T^ + -My) In - + —^ + T!^{y) + 0{e) 



6e2 3 2e 
1 1 



12e2 4e 



h{y)-^-^+T'2{y) + 0{e) 



l/2(y)ln^+.F2^(y) + 0(e) 



where T2, ^2 ^^d JF| are 0(e°) functions of y which we do not need to determine, 



A, 



Sd 



2{2n) 



[r(2 - e/2)]2r(3 - e/2)r(l + 3e/2), 



and 



/i(y) = / dv{l-v) \n{l + iyv) 
Jo 

nl 



h{y)^iy / dv 



l-v 



l + iyv 



It is straightforward to show that 



Re[/i(y)] = -5i(2;) 

^e[j2{y)] = -2S,{y)~ydyS,{y), 



(4.34) 
(4.35) 
(4.36) 

(4.37) 

(4.38) 
(4.39) 



(4.40) 
(4.41) 



where Si was defined in ( |3.9| ). We use this result, along with Eqs. (4.25) and (4.34-4.36) to write I2 (4.21) as a product 
of u^ and a series in e. In particular, we have 



e „-;2 



2/2 (w) = ^CTK ' U 



rc-2 c-2, c-2„,, (c+l)..„/^ 0.787 „„^, 

^-^ Inx + 5i y + ^ -ydySiiy) + 2.36 Ina: - 

3e"^ e e e e 



^i£^(lnx)^ 



-3(c-2)5i(y) ln2;-3(c+ 1)2/3^51(2/) Ina; + I?(2/) + 0(e) , (4.42) 



where c is given in Eq. (1.6) and 

V(y) = 0.233- (c- 2)5i(y) - {c+\)ydySM + 12 Re {^{y) + T%) ^ Tl[y)\. 



(4.43) 



V. CONNECTED PIECE OF THE CONDUCTIVITY 



The topologically distinct diagrams resulting from the expansion of the connected four-point order-parameter av- 
erage (2.18) to Oiv}) are shown in Fig. 4. Self-energy corrections are included in these diagrams since we have 
renormalized the theory, following dimensional regularization. The algebraic expressions for each allowed permuta- 
tion of wavevector and frequency in these diagrams is inserted in a^^> (2.17), thereby giving a contribution to the 
conductivity. The 0(u) diagram in Fig. 4a does not contribute to the conductivity since in this case the integral 
(2.17) separates into a product of odd integrals over ki and k2. The remaining diagrams in Fig. 4 are 0(u^)^ and 
produce 



a(4)(L^) 



128e2 
dT 



K-^ {uK-y x-'^ Re [4/b(y) +/i^^(y) +/^^^(y)] 



r(2)^ 



(5.1) 



when inserted into Eq. (2.17). The diagram in Fig. 4b is responsible for the contribution 



Ib{y) = I ki-k2 Go(ki,wi)Co(ki,wi - 22/)Co(k2,W2)C'o(k2,W2 + 22/)G'o(k3,W3)Co(ki +k2 +k3,cji +CJ2 -CJ3) 

(5.2) 



in (5.1) with y defined in (4.14) and Go and Co given by (4.22) and ( 4.24| ), respectively. The diagram in Fig. 4c 
produces the other two integrals, 



/^i)(y) ^ ki-k2 Go(ki,wi)Go(ki,2y-a;i)C'o(k2,c^2)C'o(k2,tJ2 + 2y)C'o(k3,c^3)C'o(ki+k2+k3,t^i+cj2+tJ3) 

(5.3) 



and 

-^c^^y) = / ki •k2 Go(ki,cJi)Go(ki,a;i - 2y)Go(k2, W2)Go(k2, -2y - cj2)Go(k3, a;3)Go(ki + k2 +k3,tJi +W2+W3), 

(5.4) 



(|5J) 






(a) 



(b) 



(c) 



FIG. 4. The topologically distinct diagrams in the expansion of the four-point order parameter average (2. IS) to 0{u ). 
The diagrammatic symbols are the same ones used in Fig. 3. Each diagram corresponds to several possible wavevector and 
frequency assignments, which are not shown. 



As with the integral l2{y) (4.22) for the disconnected piece, we evaluate the frequency integrals in Eqs. (5.2-5.4) 
with contour integration and use the Feynman formula (4.33) to perform the wavevector integrals. Upon e-cxpanding 
the results we have 



" ^"' 96 



1 1 , 4\ 1 /3 3 4\ 1 „ , , , ,, 0.057 _ , , ^, , 
4-2ln3J^-(4-2ln-j-Rc[A(.)] + ^+^.(.) + 0(.) 

lln^ - ^Inl Re [My)] + ^-^+TP{y) + 0(e 



0.618 



+ TP{y) + Oie) 



(5.5) 
(5.6) 
(5.7) 



where Ad and /i were defined in (4.37) and ( 4.38| ) respectively. As before, J^b, ^c and Tc are 0(e°) functions of 
y that do not need to be determined. Equations (|5.5|-5.7) are substituted into Eq. (5.1) and the result, expressed in 
terms of a product of u^ and a series in e, is 
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'^'^'\^) = ^-^>^-'u' 



2 2 2 n 088 

— r - - Inx + -Si{y) + 0.258 Inx + iilnxf - QSiiy) In x + Ciy) + 0(e) 

3e^ e e e 



, (5.8) 



where 



C{y) = 0.787- 25i(y) + - Re [4J^b(y) + J^(i)(y) + .Ff )(y)]. 



(5.9) 



In Eq. (5.8) we have used the relation (4.40) between /i and 5*1, Eq. (3.9). 



VI. ADDITIVE RENORMALIZATION OF THE CONDUCTIVITY 



The real part of th e co nductivity ( 2.15 ) is sum of the disconnected contributions, Eqs. ( 4.18 ) and ( 1.42 ), and the 
connected piece, Eq. (5.8): 



e- _,,, _ . 2n I 1 



a'{uj) = —aK.~'{l-2.Qu 



1.4 



3e2 



— In; 



e + 2cu^ 



(Inx) 



-(Inx) 



1 - (e + 2cfl^)lna; 



+-(lnx)2J 5i(y) + (e + 2cw2)(l - eln2;)52(y) + (c + 1)(77 - 2v?)ydySi{y) Inx + ^^.^(y) + ©(u^e, e^) 



where 



(6.1) 



.F(y) = 2.15i(y) - icS^iy) + -^SM + V{y) + C{y) - (c + l)ydyS2{y) 



(6.2) 



is an 0(e°) function of y. Even after renormalizing the bare quantities in the theory some poles in e remain in Eq. 
(S.l). In fact, this problem arises even in the Gaussian theory [the 1/e term in (6.1)] and indicates that we must be 
more careful when we set the cutoff A to infinity. We should write the conductivity for d < 4 as 



with 



cr'(w; d, A) = cr'(cj; d, cxi) - A{uj; d, A), 
A{u;] d, A) = a'{uj, d, cx)) - t7'(w; d, A). 



(6.3) 



(6.4) 



The cr' {lo; d, cxd) term in Eq. ( |6.3[ ) is just Eq. ( |6.1| ). By subtracting A (6.4) from (j'{uj; d, cxd) we render the conductivity 
finite in four dimensions, since we recover the theory with finite A. At low frequencies, near Tc, we expect to be able 
to approximate A by its value at Tc and to — 0: near criticality, th e IR singularities, which appear in a'{uj]d,oo), are 
absent in A since only UV physics contributes to the difference in (6^). In the minimal subtraction schem e th e poles 
of (j'^LO = 0;(i, cxd) which contain no singular temperature dependence are simply subtracted from Eq. (O). This 
situation is reminisc ent of the additive renormalization of the specific heat in the static theory ||32| . 
Inspection of Eq. (6J) gives 



A= — CTK-'(l-2.6u^ 

2r ^ 



2 
3? 



1.4_ 
e 



cu^ + — u^ + 0{%?) 



(6.5) 



and thus we write the fully renormalized conductivity cr^(w) — (t'(w) — A as 



^flH = ^'^«"^(l-2.6u2)<{-ln: 



e + 2cu2„ ,T e^ 



(lna;)2--(lna;)3 + 
6 



1 - (e + 2cv?) In a; H (In xf 



Siiy) 



+ (e + 2cu^){l - e In a;) 5*2(2/) + (c + l)(r/ - 2u'^)ydySi{y) Inx + u^Tiy) + ©(u^e, e=^) }. (6.6) 



Now we have a theory that is UV convergent as e ^ 0, XY-m odel fixed-point value u* (|2.35| ), with 77 = 2(u*)^ 
but has IR divergences as T ^ Tc and x ^ 0. Near four ( t2.36|) and, after re-summing, the series in e takes the 
dimensions, the coupling constant u flows in the IR to its form 
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<^r{^) 



2r 



-""Siiy) 



+px-PS2iy)+p'Hy) + 0{p') , (6.7) 



where p is 0(e) and is defined as 



p = e + crj 

^ 2 - d + z + 0{e^) 
2 



-{2-d + z) + Oie^). 



The (1 — 2.6m^) factor in Eq. (S.7) has been absorbed by 

changing the normahzation of (t'j^{u! = 0)^ As T ^ T^ 

and X —> terms proportional x~p in Eq. ( |6.7|) dominate 
the conductivity |3J]. From Eq. ( 4.13| ) we have 



X 



-p 



i^^r[l + Oie% 



(6.8) 



and as T ^ Tc we write 



O'fl(^) 



2r 



-aK 



P 



[l+pSi{y)+p'S2{y) + 0{p% 

(6.9) 



In Eq. (6^), the series in p coincides, to 0{p ), wit h th e 



e-expansion for the Gaussian scaUng for m 5*^ , Eq. (3 
Thus, by re-summing the series in Eq. ( |6.9[ ) we obtain, 
correct to O(e^), the Gaussian scahng form 5^, Eq. (^.6[) , 
now as a function of the critical scaled frequency y (4.14) 
and with occurrences of e replaced by p. 

The final result for complex a.c. conductivity in the 
critical regime is then (dropping the R suffix) 



a{uj) 



2r 



-an 



m 



2-d+z 



(2-d + z) 



[S{y) + 0{e^% (6.10) 



with the scaled frequency y given by Eq. (4.14) and the 
univ ersal complex scaling function S{y) given by Eq. 



VII. COMPARISON WITH EXPERIMENT 

It is instructive to compare the universal function S'(y), 
Eq. (1.7), for the critical theory (extrapolated to d~ 3), 
with both the prediction of the Gaussian theory, Eq. 
(3.7), and the experimental results of Ref. |18[. Strictly 
speaking, it is inconsistent to compare scaled data from 
different theories and experiments if the axes have been 
scaled using different exponents. However, for the sake of 
comparison, we take the viewpoint that the theory and 
experiment each determine a particular universal func- 
tional dependence S{y) and ignore exactly how S{y) and 
y are achieved. 



In this spirit, the magnitude of S{y) as a function of 
y is plotted on a log-log scale in Fig. 1 for the critical 
and Gaussian theories. Since z ^ 2 in the critical the- 
ory, the power-law behaviour at large y [a consequence of 
(1.2)] for the critical theory lies only slightly below the 
Gaussian theory. 




i<i^iy) 



FIG. 5. Comparison between the scaled a.c. conductiv- 
ity data from Booth ei al., Ref. [g_8|, on YBCO and the re- 
laxational 3D XY critical theory, (a) The scaling function 
S{y), Eq. JLtI) , using the relaxational 3D XY value z = 2.015 
(dashed curve) and using the experimental value z — 2.65 
(dotted curve) are compared with the experimental results 
(solid curves) . The magnitude of S{y) is plotted against y on 
a log — log scale. The theory is fit to the experiment using 
horizontal and vertical offsets (the horizontal offset depends 
on the value of z used), (b) The normalized phase, 2(j>{y)/n, of 
the conductivity is plotted against logj^Q y for the relaxational 
3D XY critical theory with z = 2.015 (dashed curve), the 
theory using the experimental value z = 2.65 (dotted curve) 
and experiment (solid curves). The horizontal offsets are the 
same as in (a). 

In Fig. 5, the critical theory is compared with mea- 
surements of the microwave conductivity of a thin film 
sample of YBCO in the range 45 MHz - 45 GHz near 
Tc ilq]. In this experiment, the exponent z — 2.65 ± 0.3 
and the transition temperature Tc = 89.1 ± 0.1 K, were 



determined from the power-law behaviour (1.2) expected 
at Tc- The best scaling collapse of the data determined 
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the value v = 1.0 ± 0.2 for the static exponent. In Fig. 
5a the magnitude of S{y) is again plotted as a function 
of J/ on a log-log scale. The Gaussian theory is not plot- 
ted since it lies so close to the critical theory. Since F, k 
and the prefactor of ^, which appear in both the scaled 



frequency y (4.14) and the prefactor to the conductiv- 
ity ( 6.10| ), are parameters in the TDGL theory, there is 
freedom to choose the horizontal and vertical positioning 
of the theory so as to give the best fit to the data. As 
with the Gaussian theory, the critical theory fits the ex- 
perimental scaling curve well over almost four decades in 
scaled frequency y, but deviates from the experimental 
data taken nearest to Tc. 

The dynamic exponent for the relaxational 3D XY- 
model is known to have the value z ~ 2.015. Neverthe- 
less, it is instructive to consider the z appearing in S{y), 



Eq. (L7), as an adjustable parameter. By choosing the 
experimental value z — 2.65 and adjusting the horizontal 
offset of the theory in Fig. 5, a better fit to the experi- 
mental data closest to Tc is achieved — at the expense of 
worse agreement with the rest of the data. This com- 
parison emphasizes that the experimental value z = 2.65 
seems to originate in the data set taken closest to Tc- 

The phase (j){y), Eq. (pT^), of the conductivity is plot- 
ted against logj^g y in Fig. 5b for the critical theory 
{z = 2.015), "pseudo" -theory (z = 2.65) and the ex- 
periment ilS]. As with the Gaussian theory, the critical 
theory predicts a smaller phase near Tc than seen exper- 
imentally. The "pseudo" -theory is in better agreement 
with experiment near Tc than the critical theory, but 
again does a poorer job fitting the rest of the curve. 



VIII. CONCLUSIONS 

We have examined a theory for the a.c. conductivity 
of a superconductor that includes the strong, interact- 
ing order-parameter fluctuations expecte d ne ar critical- 
ity. The FFH scaling hypothesis, Eq. (1.1), is shown 
to hold at O(e^) in the e-expansion for relaxational 3D 
XY-model critical dynamic s. T he universal scaling func- 
tion S{y) appearing in Eq. (1.1) is explicitly calculated to 
0{e^) for this dynamics, with the result given in Eq. (1.7). 
The frequency and phase behaviour expected at Tc^ Eqs. 
(1.2) and ( |l.4| ) respectively, is demonstrated. The critical 
scaling function S{y) generalizes the Gaussian result, Eq. 
(3.7), and reduces to it when z = 2. These results are 
quite general and hold, in the critical regime, for any bulk 
superconductor described by a complex order-parameter 
with relaxational dynamics. 



Since z « 2 for this dynamics, the scaling curve S{y) 
is, for practical purposes, indistinguishable from the pre- 
diction of the Gaussian theory (see Fig. 1). Therefore, in 
a measurement of the a.c. conductivity, the only indica- 
tion of a crossover from the Gaussian to critical fluctua- 
tion regime would be a crossover in the static exponent 
V. This may explain why the Gaussian theory fits the 
experimental data of Booth et al, Ref. iQ, so well over 
much of the curve in Fig. 5, even though the experiment 
is supposedly accessing the critical regime. 

The inclusion of critical order-parameter fluctuations 
in the framework of relaxational dynamics does not seem 
sufRcient to explain the deviation between the Gaussian 
scaling form and experiment |18|] observed near Tc (see 
Fig. 5). As highlighted by the fit of the "pseudo" -theory 
in Fig. 5, this deviation is connected to the large value 
z = 2.65 obtained in the experiment, which cannot ex- 
plained within any present theory [|35| . It is possible 
that this discrepancy may be due to the strong influence 
that uncertainties in the experimental determination of 
Tc have on the scaling of the data closest to Tc- More 
a.c. conductivity measurements with higher temperature 
resolution near Tc may resolve this issue, allow a more 
accurate determination of z, and provide a check on the 
scaling collapse for large y. It is also possible that the 
films studied contain strong disorder, which could affect 
the scaling near Tc- 

In this paper we have identified and dealt with the 
technical challenges involved in the organization and 
renormalization of the theory for the a.c. conductivity in 
the critical region. This work serves as a basis for examin- 
ing more complicated models, such as model F of Hohen- 
berg and Halperin pl| involving reversible couplings to a 
conserved energy-mass density field, as in superfiuid *He. 
In three dimensions z = 3/2 for model F ||30| , |36| which, 
although not observed in the a.c. conductivity dat a [18|, 
is seen in some d.c. conductivity experiments It4|, |l2|l4[ 
and simulations ^j. Another extension of the present 
theory is to consider a non-zero magnetic field, with the 
aim of examining the crossover from the zero-field critical 
scaling of the 3D XY-model to the lowcst-Landau-level 
scaling which obtains in high fields il3,pq] . 
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APPENDIX A: 



To illustrate the calculation of wavcvcctor integrals, we use this appendix to provide the details of the e-expansion 
of /|(y), Eq. (4.26|), which is reproduced, in the notation of Sec. IV E|, as 



i2iy) 



dv [l-v) / fcf 

Jkik2k3 



2 1 

af (fli + iyv)^ 



a2a3(a2 + as + 04) (as + 2iy) ' 



(Al) 



We parameterize the wavevector factors in the denominator of Eq. ( |Al[ ) in pairs using the Feynman parameterization 
(4.33), beginning with factors on the right containing ka. The denominator of the ka integral is thereby transformed 
into a quadratic form in ka and the integral is solved. The process is repeated for the remaining two wavevector 
integrals, producing 



^2(y) 



Ad 

3e 



dv (l-v) [2Jiv^0,y)-J{v,y)] 



with 



™l _-3e/2 

J{v,y)= / duidU2du^dU4:U2{l + U2Y'^uJ'^{l~UifuX'^ ^_° 3_e/2 ' 

52 9i 







where we have defined 



go = (1 - 'a4)(l + iyv) + Uigo, 

gi = 1 + U4{gi - 1), 

go = l-U3 + U3(l + 'a2){l + U2[2 + ui(l + 2^?;)]}, 
U2U3 



91 



52 



-{52[1 + Ul(l + U2)] - U2U3}, 



52 = 1 - W3 + ^2^3(2 + U2). 



In the e-expansion J in Eq. (A3) is 0(e ^) at leading order. The singularity in e is isolated by writing J as 

J{v,y) = Ja{v,y) + Jb{v,y) 
where 

Ja{v,y)^{l + iyv)-^''^^ / du2du3dui U2il + U2Y~^ ^.^u"-^ 



,^/2 



Jb{v,y) = I duidu2du3dUi U2{1 + U2Y ^ 2-^e/2 "'4 

52 



e-1 



2-£/2 "4 ' 
52 

5i 



In the e-expansion, Eq. (AlO) becomes 



1 



Ja{v,y) = - du2du3 



U2 



- ln(l + iyv) I du2du3 
Jq 



e .In gj{l + U2) 2 io gl{l + U2) 



U2 



du2du3 



U2 



.91(1 + W2) 



ln(l + U2) + - lnu3 + - ln.g2 



where J^a{v,y) is a function of w and y. The integrals in Eq. (A12) are evaluated to produce 



J„(w, y) = - In ^ - ^ In ^ X ln(l + iyv) - 0.087 + eJ^aiv, y) + ©(e^). 



The non-singular integral J;,, Eq. (All), has the expansion 

U2 



Jb{v, y) = duidu2du3dui 



1 



gl{l + U2)ui 



~9l 



- 1 



+ eTb{v,y)+0{e') 



0.103 + e.Ffc(«,2;) + 0(e2), 



(A2) 

(A3) 

(A4) 
(A5) 
(A6) 

(A7) 
(A8) 

(A9) 

(AlO) 
(All) 



+ eTa{v,y) + 0{e^), (A12) 



(A13) 



(A14) 
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where J^b{v,y) is a function of v and y. By combining the Eqs. (A13) and (A14) in Eq. (A9) we may use this result 
for J in If, Eq. (|A^), to obtain the result quoted in Eq. (4.34), with 



•^2 (y) = I 



dv{l - v){2[Ta{0, y) + ^6(0, y)] - [Ta{v, y) + Tb{v, y)]}. 



(A15) 
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